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Abstract 

Bottom quark-induced processes are responsible for a large fraction of the LHC discovery 
potential, in particular for supersymmetric Higgs bosons. Recently, the discrepancy between 
exclusive and inclusive Higgs boson production rates has been linked to the choice of an 
appropriate bottom factorization scale. We investigate the process kinematics at hadron col- 
liders and show that it leads to a considerable decrease in the bottom factorization scale. 
This effect is the missing piece needed to understand the corresponding higher order results. 
Our results hold generally for charged and for neutral Higgs boson production at the LHC 
as well as at the Tevatron. The situation is different for single top quark production, where 
we find no sizeable suppression of the factorization scale. Turning the argument around, we 
can specify how large the collinear logarithms are, which can be resummed using the bottom 
parton picture. 



L HIGGS BOSONS AT THE LHC 

The combined LEP precision measurements [1] suggest the existence of a light Higgs boson. In the 
case of a single Standard Model Higgs boson the LHC promises multiple coverage for any Higgs boson 
mass, which will enable us to measure its different decay modes and extract the couplings [2]. For a 
supersymmetric Higgs sector this coverage has to rely on fewer Higgs boson decay channels [2,3]. This 
is a direct consequence of the structure of the Higgs sector: while the Minimal Supersymmetric Standard 
Model (MSSM) predicts a light Higgs boson, it also predicts an enhancement of the coupling to down-type 
fermions, at the expense of the branching fractions to gauge bosons. This enhancement is an outcome from 
the two Higgs doublet structure in the MSSM: one doublet is needed to give mass to up-type, the other to 
down-type fermions. The vacuum expectation values of the two doublets are different, parameterized by 
tan/? = V2/V1. In addition to a light scalar Higgs boson, the two Higgs doublet model includes a heavy 
scalar, a pseudoscalar, and a charged Higgs boson. None of these additional particles have a mass bounded 
from above, apart from triviality or unitarity bounds. 

Of course, observables linked to properties of a light Higgs boson can serve as a probe if a new scalar 
particle is indeed consistent with the Standard Model Higgs boson [4,5]. There is, however, only one way 
to conclusively tell the supersymmetric Higgs sector from its Standard Model counterpart: to discover the 
additional heavy Higgs bosons and determine their properties. 
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At the LHC, the possible enhancement of down-type fermion Yukawa couplings by powers of tan (3 can 
render the search for a heavy scalar and pseudoscalar Higgs boson promising. For small values of tan (3 
the Yukawa coupling of the charged Higgs boson is governed by the top quark mass m t j tan/5, whereas 
for larger values of tan (3 the bottom Yukawa coupling m b tan (3 dominates. While the chances of finding a 
heavy Higgs boson with a small value of tan (3 at the LHC are rather slim, the discovery of all heavy Higgs 
scalars in the large tan f3 > 10 regime is likely. This reflects the fact that the reach of the LHC for charged 
and neutral Higgs bosons is to a large degree owed to scattering processes which involve incoming bottom 
quarks. The completely exclusive processes are 

99 -> btH~ gg -> 66$, (1) 

where $ = h°, H°, A denotes any neutral scalar Higgs boson. From an experimental point of view, the 
identification of these final state bottom jets is tedious, because the dominant contribution comes from phase 
space configurations where the incoming gluons split into two collinear bottom quarks. The bottom quark 
rapidity distribution peaks at rapidities around two and the transverse momentum distributions around the 
bottom quark mass. After adding in the efficiency for a bottom tag this becomes a heavy price to pay in the 
analysis. Therefore one usually prefers to look for more or less inclusive channels 

gb^tH- gb^b<$> 66^$. (2) 

At this point we do not explicitly discuss the bottom-induced inclusive channels bb — > W + H~ [6] and 
bg — > tW~ [7], which are both known to next-to-leading order QCD, because their impact on the discovery 
potential of the LHC is not drastic. We emphasize, however, that our argument will hold for them the same 
way we apply it to the processes in eq.(2). 

All bottom-inclusive channels suffer from an additional uncertainty: the choice of the factorization scale 
of the bottom parton. In contrast to the gluon density, the dependence of the bottom parton density on the 
factorization scale is large even for scales above (9(100 GeV). Recently, it has been observed in higher 
order calculations [6-10], that varying the bottom factorization scale around a smaller central value yields 
a more stable perturbative behavior. The same choice of scales can resolve the discrepancy between the 
total inclusive and exclusive rate, which is most prominent for the production channel bb — > $. However, 
choosing scales according to perturbative behavior is difficult, because most processes include inherent 
cancellations between different contributions, and it is hard to define which contributions actually have to 
be stable. Instead of this somewhat soft argument, we will derive an appropriate bottom factorization scale 
from the kinematics of the exclusive production process. In the following two sections we first investigate 
charged Higgs boson production, because it involves only one incoming bottom quark and two heavy central 
decay products. In Section V we will then generalize our result to neutral Higgs boson production and 
compare it to single top quark production. 

Conventions: Throughout this paper we show consistent leading order cross section predictions, includ- 
ing the respective one loop coupling constant, running heavy quark masses, and CTEQ parton densities. 
Unless stated otherwise, we assume tan/3 = 30 for all MSSM processes. The exclusive cross sections are 
quoted with a massive 4.6 GeV bottom quark in the matrix element and the phase space, while the bottom 
Yukawa coupling is set to the running bottom quark mass. 
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II. BOTTOM PARTON SCATTERING AT THE LHC 



Heavy flavor-induced search channels for supersymmetric Higgs bosons have been explored for many 
years [11]. To begin with the charged Higgs boson, three search modes have been investigated: (1) charged 
Higgs bosons can be pair produced in a Drell-Yan type process, mediated by a weak interaction vertex [12]. 
Moreover, they can be pair produced at tree level in bottom quark scattering [6] or through a one loop 
amplitude in gluon fusion [13]. (2) One charged Higgs boson can be produced together with a W boson via 
scattering of two bottom quarks or in gluon fusion [14]. (3) The charged Higgs boson can be produced in 
association with a top quark, which seems to be the most promising search channel [7,8,15,16]. The charged 
Higgs boson can be detected either decaying to a top and a bottom quark [17] or decaying to a tau lepton and 
a neutrino [18]. Both LHC collaborations have tried to reproduce these phenomenological analyses, most 
successfully in the case of the decay to tau leptons [19,20]. However, these analyses have to be taken with 
a grain of salt, if they rely on the standard Monte Carlo generators, because some do not include a running 
bottom Yukawa coupling. From the next-to-leading order QCD calculations [7,8] we know that using the 
bottom pole mass as the Yukawa coupling severely overestimates the rates, as one would expect from what 
we know about Higgs boson decays to bottom quarks 1 . 

The case of neutral scalar Higgs bosons in the MSSM has been discussed in similar detail. In the Stan- 
dard Model the inclusive production process bb — > $ is a small correction to the inclusive gluon fusion 
channel [21]. In a supersymmetry context, and in particular for large tan (3, the bottom quark-induced 
process dominates gluon fusion [2,10,22,23]. The typical Higgs decays are the same as for a light Stan- 
dard Model Higgs boson, except for heavy MSSM scalars, where decays to muon or tau pairs are most 
promising [24]. Additional problems occur in the so-called intense coupling regime, i.e. the region with an 
intermediate pseudoscalar Higgs boson mass and strong mixing effects [3,25]. The two scalar Higgs bosons 
can be detected and possibly be resolved in weak boson fusion with a subsequent decay to r pairs [3]; 
unfortunately, for large values of tan (3 the mass splitting becomes too small to distinguish the two mass 
peaks. In that region, the most promising way to search for and separate the heavy Higgs states is the pro- 
duction process involving b jets with the Higgs bosons decaying to muons [26]. To distinguish the gluon 
fusion process from the bottom quark-induced production, one can require observation of only one final 
state bottom quark, i.e. using the partly inclusive channel gb — > 6$ [9,27]. 

Motivated by this vast number of analyses in the MSSM Higgs sector we turn to the simplest process 
available: charged Higgs boson production in association with a top quark involves only one incoming 
bottom quark and is an appropriate starting point to understand the issue of bottom partons at the LHC. The 
features of the exclusive and the inclusive charged Higgs boson production process 

gg - UH- gb - tH~ (3) 

have been investigated in detail in Ref. [8]. Let us briefly review the most important observations: 

For reasons described above, the searches for charged Higgs bosons (decaying to tau leptons) at the LHC 
usually do not require the observation of a final state bottom quark. The exclusive process gg — > btH~ is 



We can compare these next-to-leading SUSY-QCD calculations [7,8] to the usual heavy parton subtraction schemes 
which add the exclusive and the inclusive channels [11]: for large final state masses, the higher order calculation of 
the inclusive channel is just the perturbatively consistent extension of the latter. 
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then dominated by collinear splitting of one of the gluons into a bottom quark pair: 
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These logarithms in the total cross section can be resummed, which is precisely the definition of bottom 
partons [11]. The factorization scale of the bottom partons is defined as the maximum transverse bottom 
quark momentum up to which the asymptotic form of the cross section is assumed to hold, and up to 
which the logarithms log p T , b /m b are then resummed. This means that in eq.(4) one can identify p™^? 
with ji Fb . If we assume p™^ x ^> m b , which as we will see is true for all processes we consider, the 
numerical value for the bottom quark mass will have no impact on our argument — we could neglect 
it altogether [8,9,11]. The factorization scale is per se an artificial parameter which has to vanish after 
including all orders of perturbation theory. In the case of neutral Higgs boson production bb — > $ this has 
recently been demonstrated, including the NNLO QCD corrections [10]. We are lucky in the case of bottom 
quarks: the comparison between the actual and the asymptotic form of the exclusive cross section allows us 
to estimate and therefore the bottom parton factorization scale for the inclusive process. This also tells 
us how large the 'large logarithms', which are resummed using bottom parton densities, actually are. 

There are two caveats to be kept in mind, though. First, the inclusive process gb — > tH~ assumes that 
the entire cross section comes from a phase space region in which we can neglect the transverse momentum 
of the bottom jet appearing in the exclusive process gg — > btH~ . This is probably a good approximation, in 
particular after including detector effects, but one has to be aware of it. Second, naive dimensional analysis 
suggests np,b = mt + m H, which is often used. This does not have to be correct. The only thing dimensional 
analysis tells us is /j, Ftb oc (m t +m H ), as long as the production process is dominated by the threshold region. 
It has been shown that the proportionality factor does not at all have to be unity and that a wrong choice of 
scale leads to a systematic overestimate of the cross section, as it is obvious from eq.(4) [8,9,28]. 



III. ASYMPTOTIC BEHAVIOR IN THE BOTTOM QUARK VIRTUALITY 

To make our argument as transparent as possible, we will investigate the relation between the bottom 
factorization scale and the threshold mass in two steps. The threshold mass we will refer to as M; for the 
charged Higgs boson production process this means M = m t + m H . We can generally rewrite the exclusive 
and the inclusive production processes including one bottom parton as 

gg ->■ bX M bg ->■ X M , (5) 

where Xm denotes the heavy final state particles. For a typical gluon-induced LHC production process 
close to threshold we expect the invariant mass of the heavy system and the threshold mass M to be similar. 
In the case of gg — > btH~ we investigate how close to threshold the production takes place in Fig. 1. The 
parameter AM is defined as the difference between the invariant mass of the heavy tH~ system and the 
threshold mass M = m t + m H . We indeed see that the distributions peak at small values AM/M < 1/8, 
even though there are sizeable tails towards larger invariant masses. 

In the first step of our argument, we investigate the maximum value for the intermediate bottom quark 
virtuality <5™ ax , up to which the asymptotic form of the exclusive cross section holds 

Q™ ax = C Q M. (6) 



4 




200 400 2 4 

AM [GeV] Pz.b^T.b 



Figure 1. Normalized distributions for the hadronic charged Higgs boson production gg — > btH~. Left: dif- 
ference between the invariant mass of the tH~ system and its threshold mass M = mt + mn- Right: ratio of the 
longitudinal and transverse momentum of the bottom jet p = p z ,b/PT,b in the partem center-of-mass system. The 
steeper set of curves is after a cut Q b > M/5. 

In Section IV, we then estimate how far the asymptotic form in terms of the bottom transverse momen- 
tum is valid. From Ref. [8] we expect a relation p™ b x ~ M/6. If we assume that this reduction will consist 
of the suppression from the asymptotic behavior in the virtuality and an additional suppression when we 
move to the asymptotic behavior in the transverse momentum, we obtain 

lM Fjb = p^f = C p Qr x = C P C Q M (7) 



To understand the asymptotic behavior of the hadronic cross section a(pp — > J)Xm) as a function of 
the bottom quark virtuality, we rewrite the integration over the phase space and the parton momentum 
fractions [9]. As long as we are interested in the behavior of the cross section for large values of the rapidity 
we can safely neglect the bottom quark mass: 

l l r s - M2 o r s /-! lo §# l 

,7 = —-/ dQ 2 b ds I dyC gg - 2 \M\* (8) 

Here Q b is the intermediate bottom quark virtuality, y the rapidity, and y/S and ^/s are the hadronic 
and partonic collider energies. The factor 1/s 2 in the integrand is obvious from the difference in mass 
units between the matrix element and the partonic differential cross section da/dQl, as it originally ap- 
pears in the integral. The parton densities are denoted as C = P i / p (xi)Pj/ p (x2). At this point we make 
a few simplifying approximations: since we want to show that the asymptotic behavior in the virtuality 
is a process-independent phase space effect, we neglect all structures in the matrix element, except for 
the asymptotic behavior in the virtuality. The asymptotic form of the differential hadronic cross section 
da/dQb oc 1/Qb translates into \A4\ 2 = S 2 a /Ql. The factor S 2 we introduce to absorb units of energy, it 
could as well be M 4 . For reasons which will be obvious later, it could not be a partonic variable, since we 
need to keep track of the powers of the parton momentum fraction. Furthermore, we assume that the steep 
gluon densities balance each other for the production of heavy states Xi = x 2 = \fx for x 1: x 2 1. The 
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Figure 2. Ratio of partem distributions in the proton and the function x~ n for two different values of the factor- 
ization scale. All curves are normalized to their values at x = 0.1. The different lines correspond to the gluon (solid, 
n = 2), down-quark (dashed, n = 1.1), anti-up-quark (dotted, n = 1.7), and bottom quark(dash-dotted, n = 2) 
content. We use CTEQ6L partem densities [29]. 



approximate hadronic cross section now reads 



2a [S-m* dQ b yi dx 
2a C fS~ M2 dQ b 



16n Jo Q b 



mot)) (9) 



using 

p u = - J! % why [ 1 -^ + < 1 -^-' to H~ *ss- do) 

with x = s/S and r = (Ql + M 2 )/S. The function F(t) is as a correction to the asymptotic behavior 
of the virtuality distribution da/dQ b oc 1/Q b . In an intermediate step we have approximated the incoming 
parton luminosity by a simple power suppression C = C /x j ~ 2 . As a general parameterization of the 
partem densities, this is certainly not a good idea. However, looking at the production of heavy particles 
at the LHC we probe momentum fractions between 1CT 1 and few times 1CT 3 . In Fig. 2 we show different 
parton distributions multiplied by x n for two values of the factorization scale. The ratio is normalized to its 
value at re = 0.1. Looking at the different values of n we see that our argument will at this point become 
dependent on the parton the incoming bottom parton sees on the other side. On the other hand, we also see 
that for gg and gb initial states the approximation C oc l/(xiX2) 2 works very well. From Fig. 2 we obtain 
j — 4 for gluon or bottom initial states in the definition of F(t). In Fig. 3 we show the behavior of F(t) as 
a function of the bottom quark virtuality, the way it looks for a 500 GeV charged Higgs boson at the LHC. 
The case of j ; = 2 corresponds to constant parton densities C = Cq. 

Let us now turn to a detailed discussion of the function F(t), shown in Fig. 3. As mentioned above, 
F(t) is a correction to the known asymptotic behavior of the differential hadronic cross section with respect 
to the bottom quark virtuality. As in Ref. [8], we show the (normalized) curves for Q b da/dQ b , which we 
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Figure 3. The normalized function F(t), defined in eq.(10). The hadronic center-of-mass energy is set to 
f S = 14 TeV and the threshold mass to 675 GeV, corresponding to a 500 GeV charged Higgs boson. We dis- 
play the behavior of the plateau in Q b for different values of j, which arise from the x±x 2 behavior of the partonic 
cross section. 

numerically obtain for the exclusive btH~ production process, in Fig. 4. The first thing we notice from the 
exact results in Fig. 4 is that, as a function of Qb/M, the curves for different Higgs boson masses are almost 
identical. The only major difference arises from the finite bottom mass effects, since their onset does not 
scale with Q b / M, but with m b /Q b . On the other hand, these mass effects are understood and of no relevance 
to our argument, which is concerned with the upper end of the asymptotic behavior in Q b . We roughly 
indicate this upper end of the plateau with a dotted line, where Q b da/dQ b has dropped to half of its plateau 
value. In the first row in Fig. 4 we observe how the plateau in the virtuality does not at all extend to Q b = M. 
The asymptotic approximation of the virtuality distribution is valid only up to values Ql 1 ^ ~ M/2.5. In the 
second row we see that this picture changes when we ignore the gluon densities, but keep everything else, 
like in the complete numerical analysis: now the asymptotic behavior extends to Q b > M. In other words, 
the short plateau in Q b da/dQ b is an effect of the steep gluon density in the proton. In the last row we also 
show how the approximation of the gluon luminosity C gg {x\, x 2 ) ~ 1/ (xix 2 ) 2 , which is a major ingredient 
we use to derive the approximate form F(t), works very well for the virtuality distribution. 

Since our main interest is the size of the bottom quark virtuality plateau, we need to compare the ap- 
proximate form of F(t) in Fig. 3 with the exact distribution Q b da/dQ b in the first row of Fig. 4. In all 
figures we include the line which indicates where Q b da/dQ b has decreased to half the plateau value. This 
is as good a measure for the extension of the plateau as any other. For the discussion of the approximate 
function F(t) we prefer a better suited measure: the turning point of F(r (log Q b /M)) 



d 2 F{r{QD) 



d(\ogQ b y 



= 0. (11) 



In Fig. 3 and Fig. 4 we see that this definition gives essentially the same Q™ ax values as the dotted line. The 
numerical values for Q™ ax which we compute from the definition of F(t) are given in Tab. I for different 
values of j. The case j = 2 corresponds to the case with constant parton densities. It is in very good 
agreement with what we see in the second row of Fig. 4. The case j = 4 should give the extension of 
the plateau for gg- and gfr-initiated processes, for example the charged Higgs boson production process. 
We see that the approximation \M\ 2 = S 2 a /Ql does not give a perfect prediction of Q b iax , as it leads 
to Qf iax ~ M/1.8. For interfering s and t channel diagrams in the production process gg — > bX M , the 
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LHC: gg^btH" m H =250,500,1000 GeV 




LHC: gg->btH' with P s (x)=l 




LHC: gg->btH" with P (x)=x 




1 Q b /(m t +m H ) 10 1 p Tb /(m t +m H ) 



Figure 4. Normalized distributions for the hadronic charged Higgs boson production process, for the com- 
plete gluon density (first row), for a constant gluon density (second row), and for the approximate gluon density 
P(x) = l/x 2 (third row). The left column shows the bottom quark virtuality distribution, the right column the bottom 
quark transverse momentum. The normalization for the largest Higgs boson mass is by the total rate; for all other 
masses the curves are normalized such that their maxima coincide. The normalization factors for the virtuality and 
the transverse momentum are identical. We note that for a comparison with the approximate form F(t) we have to 
identify M = m t + mn- 

common denominator in the differential cross section typically becomes \/{sQl), while the numerator is 
dominated by the heavy mass. In our simple approximation we did not take into account this additional 
factor l/s, which increases j to 5 and moves the turning point to Q™ ax < M/2 and therefore much closer 
to the values seen in Fig. 4. 
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Vs 


M 


gmax 


Qr x /M u = 4) 


gmax 


Qr x /M (j = 5) 


gmax 


Qr x /M (j = 2) 


2000 


130 


74 


1/1.76 


64 


1/2.03 


123 


1/1.06 




250 


142 


1/1.76 


123 


1/2.03 


235 


1/1.06 




500 


282 


1/1.77 


245 


1/2.04 


463 


1/1.08 




1000 


559 


1/1.79 


488 


1/2.05 


908 


1/1.10 


14000 


130 


72 


1/1.81 


64 


1/2.03 


119 


1/1.09 




250 


138 


1/1.81 


121 


1/2.07 


221 


1/1.13 




500 


271 


1/1.85 


238 


1/2.10 


416 


1/1.20 



Table I. Maximum values for the bottom quark virtuality at the Tevatron and at the LHC, as defined in eq.(10) and 
eq.(ll). The values of j correspond to the power of x in eq.(10), as it arises from the x\X2 behavior of the partonic 
cross section. 

To summarize this section: we have shown that partonic phase space effects are responsible for the 
maximum value Qjf ax , up to which the asymptotic behavior of the hadronic cross section with respect to the 
bottom quark virtuality is valid. Our very simple approximation agrees with the numbers we obtain for the 
full hadronic process gg — > btH~ in Fig. 4: 
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Qr x = c Q M <v„~ 2 c q^~—- (12) 
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IV. ASYMPTOTIC BEHAVIOR IN THE BOTTOM QUARK TRANSVERSE MOMENTUM 

From the discussion of the bottom quark virtuality one would expect to be able to follow a similar set 
of arguments for the bottom quark transverse momentum. Unfortunately, the phase space parameterization 
reflects the fact that the hadronic cross section factorizes in the virtuality and not in the transverse momen- 
tum. Instead, we choose a different path: we know that the asymptotic behavior da/dpT,b oc l/pr,b nas 
to hold for small transverse momenta. In that regime the longitudinal momentum of the outgoing bottom 
quark in the center-of-mass system will be much larger than the transverse momentum: p = p Z: b/pr,b 3> 1. 
On the other hand, we want to push this approximation to its limits, pel. The general relation between 
the virtuality and the transverse momentum of the bottom jet for the production of a heavy system X M at 
threshold is: 

(.ir7p_ p \„PT* (13) 



M 2 M 2 V v / M 

In the intermediate step we have assumed that in the limit of large bottom quark virtuality and large trans- 
verse momentum the bottom jet has dominantly a transverse momentum direction, i.e. that the bottom jet is 
central in the detector. We now require that the phase space region which forms the upper end of the plateau 
in the virtuality (<5jf ax ) is also responsible for the upper end of the transverse momentum plateau (p™^ x ). 
From the approximate result in eq.(12) we obtain 

p F ,b = P^f = C p Qr x = C p C Q M C p =C Q =i (14) 

" nl "'" v approx 2i 



approx 
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To understand this effect in more detail, we now keep p as a free parameter and assume that the heavy 
system X M be produced at threshold. We find: 



s-M 2 



(15) 



PT,b = 



We can solve the equation for y/s and insert it into eq.(13), which leaves us with: 



M 2 M 2 



Ql _ Pr,b 



( 



1 + p 2 - p 




M 2 



(16) 



This function has two limiting cases: for small transverse momentum (p 3> 1), p Tjb scales with the 
virtuality p Tb ~ Q b . This means that the plateau in the transverse momentum will extend to the same value 
as the plateau in the virtuality. This observation suggests that (at least for not too large transverse momenta) 
the phase space region which dominates the high end of the plateau in Q b will also be responsible for the 
upper end of the plateau in p Tjb . We have confirmed this assumption for charged Higgs boson production 
explicitly. In contrast, for very large transverse momentum (/) < 1), a correction to the linear relation be- 
tween pT,b and Q b occurs. This correction pushes pTf, to smaller values, and the corrections becomes bigger 
for small values of pr,b/M, which is what we expect from the shifted and softened transverse momentum 
plateau in Fig. 4. 

To translate Q™ ax into p™ b x we evaluate eq.(16) numerically. The left hand side we substitute by Cq, 
which according to the previous section assumes numerical values of 1/2 to 1/3. Fig. 5 shows the cor- 
responding values of C p = PtT/QT^ f° r different p values. Indeed C p — 1 holds for small transverse 
momentum p > 5. For all other values of p the transverse momentum is always considerably smaller than 
the virtuality. In other words: if we want the transverse momentum to be as large as possible for a given 
virtuality, we would have to make p large, i.e. make the longitudinal momentum even larger. The limiting 
factor will be once more the steep gluon luminosity. In the opposite regime p <C 1, which will be preferred 
by the gluon density, we find a very substantial reduction factor C p < 1/2. 

In the right panel of Fig. 1 we plot the p distribution for exclusive charged Higgs boson production. It 
is clearly peaked at small values of p, which means large values of pr,b- The peak becomes considerably 
more pronounced if we only allow virtualities above M/5, which limits the phase space to the transition 
region in the virtuality plateau. According to Fig. 5, the region p < 1 leads to transverse momenta much 
smaller than the virtuality. Even though the actual distribution in Fig. 1 peaks at p = these events will not 
anymore contribute to the upper end of the plateau in p T:b ; instead the plateau in p T>b will be softened. At 
the opposite end, p > 2, we would expect a negligible suppression factor C p and large transverse momenta 
from virtualities around Q™ ax , but the actual distribution shows how the gluon luminosities cut heavily into 
this region. What is left is the intermediate region p = 0.5 ■ • ■ 1.5, which interpolates between the extremes 
and contributes most to the upper end of the plateau in pT,bda/dpT,b- Even though the numerical details used 
in this argument are not process independent, the fact that one has to interpolate between the two extreme 
regions in phase space is completely general. 

This leaves us with two conclusions concerning the connection between the asymptotic regions in the 
bottom quark virtuality and transverse momentum. First, the plateau will not just translate from Q b da/dQ b 
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Pz,t/PT,b 



Figure 5. The numerical solution of eq.(16) for two different values of Cq = Qb/M, giving C p = pr,b/Qb as a 
function of p = p z ^/PT,b- The longitudinal momentum of the bottom quarks is defined in the parton center-of-mass 
system. 

to PT,bda/dpT,b- Instead, it will be softened. Second, we can extract our approximate prediction for C p from 
Fig. 5 (keeping in mind Cq ~ 1/2.5) and compare it with what we find in Fig. 4: 



1 1 



^F,b = PtT = C p QT™ C p ~ T~T " ' o C p ww ~ o ( 17 ) 

approx 1.4 Z 1:111 
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Together with the results presented in Section III we now understand that for charged Higgs boson 
production the plateau in pr,bdo j dp F ,b does not extend to values p™^ x ~ M. Making use essentially of 
phase space effects we instead find p™ b x ~ M/4. This result confirms the value p™^ x ~ M/5, which we 
find directly from the exclusive process gg — > btH~ [8], which implies that using the naive bottom quark 
factorization scale fi F:b = M will overestimate the cross section considerably. On the other hand, higher 
order QCD corrections [6-10] soften the dependence on the factorization scale considerably. 

We can also turn this argument around: the bottom parton approach means integrating over the additional 
bottom quark phase space and resumming the logarithms including the transverse momentum. Going back 
to eq.(4), we see that the terms which we resum are at maximum a s \og[(p^ b x /m b ) 2 ]. For a charged Higgs 
boson of mass 500 GeV this gives 6.7a s . For a threshold mass M = 130 GeV the logarithm only yields 
moderate 3.5a s , using p™^ = M/5. 



V. SIMILAR AND NOT SO SIMILAR PROCESSES 

A. Neutral Higgs Boson Production 

In the previous sections we derived the appropriate bottom parton factorization scale for the associated 
production of a charged Higgs boson and a top quark. For different reasons this process is particularly well 
suited for the bottom parton description: there is only one bottom parton, the top-Higgs system is very 
heavy, it will be produced close to threshold, and (for those reasons) it will be slow-moving and central in 
the detector. A process which is particularly important for light as well as for heavy supersymmetric Higgs 
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Figure 6. Normalized distributions for hadronic neutral Higgs boson production gg — ► bb<& and single top quark 
production qg — ► frig'. Left: difference between the invariant mass of the heavy system and its threshold mass. Right: 
ratio of the longitudinal and transverse momentum of the bottom jet p = p z ,b/PT,b in the parton center-of-mass system. 

bosons is the production of a neutral scalar in association with two bottom quarks, which at the different 
level of inclusive versus exclusive description reads [22,10]: 



gg -> bb§ 



bg — > M> 



bb -> $ 



(18) 



A phenomenological interesting aspect is that these kind of channels with at least one tagged bottom 
jet [9,27] can prove an enhanced bottom Yukawa coupling. We know that just like in the charged Higgs 
boson case, the factorization scale of the bottom parton has to be chosen well below [8,9]. The first 
reason is that again a heavy system (X M = 6$) is produced close to threshold, Fig. 6. In Fig. 7 we show 
the plateau in Q b da/dQ b , similarly to the figures in Ref. [9], and obtain Q™ ax ~ M/2.5. Like for the 
charged Higgs boson case we see that the curves are nearly degenerate for different Higgs boson masses, 
i.e. Ql 11 ^ oc M. For the neutral Higgs bosons an additional curve is included in Fig. 7, assuming a 130 GeV 
light scalar MSSM Higgs boson. In this case the plateau is not particularly wide, since the bottom mass 
effects bend down the curves at fairly large values of Qb/M. On the other hand, these effects are understood, 
so that the altered shape can be treated just like a plateau. Moreover, the dominant effects will arise from 
the upper end of the curve, where the logarithms are largest. 

Comparing Fig. 4 with Fig. 7, we see that the extension of the asymptotic behavior in the bottom quark 
virtuality at the LHC is just the same for charged Higgs boson production (bg — > tH~) as for neutral Higgs 
boson production (bg — > 6$), for similar heavy state masses. An aspect we did not discuss in the charged 
Higgs boson case is that the plateau seems to extend to slightly larger virtualities for smaller threshold 
masses, in particular for a 130 GeV neutral Higgs boson. The reason is that, in general, Xm will be 
produced relatively closer to threshold for heavier states, i.e. AM/M becomes smaller, even though Fig. 1 
and Fig. 6 show that in absolute numbers AM becomes slightly larger for heavier states X M - 

Up to this point we did not have to specify the collider energy in our approximation. Moreover, charged 
Higgs boson production is only relevant at the LHC. The distributions for neutral Higgs boson production at 
the Tevatron are given in Fig. 8. While for small Higgs boson masses of 130 GeV the bottom quark virtuality 
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LHC: gg^bbH m H =130,250,500,1000 GeV 




Figure 7. Normalized distributions for the hadronic neutral Higgs boson production and for exclusive single top 
quark production at the LHC. The normalization for the largest Higgs boson mass is by the total rate; for all other 
masses the curves are normalized such that the virtuality distributions coincide at their maxima. The normalization 
factors for the virtuality and the transverse momentum are identical. The dashed curve uses a mathematical cutoff 
0.46 GeV for the bottom quark mass and 130 GeV for the Higgs boson mass. It is normalized to match the curve for 
the physical bottom quark mass for large virtuality. The curves are ordered on their downward slopes by decreasing 
final state mass towards larger values of Q™ ax and p™^ x . In the case of single top quark production we use three 
different top quark masses to illustrate the dependence on the final state mass, even though the top mass is a measured 
parameter . 



plateau looks very similar to corresponding curves for the LHC, our picture starts to break down for very 
large Higgs boson masses. In that case the limited hadronic collider energy does not allow production 
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Tevatron: gg^bbH m H =130,250,500,1000 GeV 




Figure 8. Normalized distributions for the hadronic neutral Higgs boson production at the Tevatron. The normal- 
ization for the largest Higgs boson mass is by the total rate; for all other masses the curves are normalized such that 
the virtuality distributions coincide at their maxima. The normalization factors for the virtuality and the transverse 
momentum are identical. The dashed curve uses a mathematical cutoff 0.46 GeV for the bottom quark mass and 
130 GeV for the Higgs boson mass. It is normalized to match the curve for the physical bottom quark mass for large 
virtuality. The curves are ordered on their downward slopes by decreasing final state mass towards larger values of 

gmax an( j p max 

of heavy states plus a bottom jet with sizeable virtuality /transverse momentum. Of course, the case of a 
500 GeV Higgs boson at the Tevatron is phenomenologically irrelevant. On the other hand, we learn that 
implicitly we have used another approximation in our discussion of the LHC processes: that the production 
rate even for large virtuality is limited only by the parton densities and the parton energy, never by the 
hadronic collider energy. Implicitly we check this requirement in Fig. 4, when we confirm <5™ ax ~ M after 
neglecting the effect of parton densities. 

Last but not least, the question is what will happen for the completely inclusive Higgs boson production 
(bb — > $)? From Ref. [9] we know that the argument for the plateau in Qbda/dQb works just the same 
way as before, except that we now use the semi-exclusive process bg — > 6$ to compute Qjf ax . Using 
the approximation described in Section III we understand what happens: the mass of the heavy system is 
still M = m$, but the gluon which splits into a bottom quark pair now sees a bottom parton on the other 
side instead of the gluon in the completely exclusive process. In Fig. 2 we see that the approximations 
£>bg ~ l/^i 2 ^) 2 and x\ ~ x-2 still work fine. In a way this is a consistency check, of course. A significant 
change in the behavior of the bottom partons between the first bottom jet and the second bottom jet to be 
integrated over in the process gg — > 66$ would pose a serious problem for factorization in general. This 
result is in very good agreement with the results of the NNLO QCD calculation for this process: the small 
bottom quark factorization scale indeed yields perturbative corrections which are well under control [10]. 

The second step is again the transition from the plateau in virtuality to that in transverse momentum. 
Not surprisingly, from the exclusive process we numerically find the same behavior as for the charged Higgs 
boson, Fig. 4. Our argument in Section IV also works just the same way for neutral Higgs bosons. From 
the exclusive process we numerically find fj, F>b ~ M/5, which we can understand in complete analogy to 
the charged Higgs boson case in the previous sections. 
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B. Single Top Quark Production 



In contrast to the Higgs boson production mechanisms, the bottom parton picture in single top quark 
production [30,31] 

ug — > btd (ub — > td) dg — > btu (db — > tu) (19) 

has never posed a problem. The difference between the two above processes is that the first one will involve 
valence quarks at the LHC, while the second one will not. Looking back, Fig. 2 shows that these channels 
should look slightly different, if our argument in Section III is correct. 

In the left column in Fig 7 we see how the single top quark case differs from the Higgs boson production. 
First of all, the plateau in the virtuality extends considerably further, typically to Q™ ax ^ M/1.5. This is in 
agreement with the less steep parton densities for the quarks, which the splitting gluon sees. If we have a 
closer look, we even see that the plateau extends further in the case of an incoming valence quark than for 
a sea quark, which is in agreement with the approximate parton densities, Fig. 2. 

However, our approximation has to be looked at with some care, since now we cannot assume x\ ~ x 2 
anymore. Furthermore, in Fig. 6 we see that the single top quark production does not at all happen close 
to threshold. It peaks around AM ~ myy, which reflects the fact that one could integrate over the phase 
space of the outgoing jet and regard the single top quark process as bW scattering. This gives the outgoing 
jet a transverse momentum kick of the order of the W boson mass. The fact that the invariant mass of the 
heavy system X M = tj, as it appears in eq.(9), is easily twice the threshold mass, again contributes to the 
larger values of Qf mx . Finally, as discussed in the context of neutral Higgs bosons, production away from 
threshold lifts the degeneracy of different values of Q™ ax /M, pushing the lighter states to higher values of 

/Omax 

Vb • 

We will not discuss the transition from the bottom quark virtuality to the transverse momentum in any 
detail. From Fig. 7 we see that the plateau is softened and p ma , x ~ Q™ ax /2, as in all other processes dis- 
cussed before: the interpolation argument presented in Section IV describes the single top quark production 
perfectly well. This part of our argument indeed holds independently for all processes considered in this 
paper. 

Recently, a similar issue of bottom partons was discussed [32] in the framework of single top quark 
production at a linear collider — > vtb. The authors find sizeable differences between the finite 
prediction and the (massless) structure function approach, predominantly close to threshold -Js < m t + 
m b + 10 GeV. These differences can in part be traced back to phase space effects. From Fig. 1 and Fig. 6 
we see that this region of phase space contributes little to the Higgs boson sample at the LHC, after we 
convolute the partonic cross section with the gluon densities, integrating over the entire partonic energy 
range. It will have even less impact on the total rate once a minimum transverse momentum of the Higgs 
boson decay products is required. While for a linear collider the bottom quark mass is an important source 
of theoretical uncertainty and the collinear logarithms (multiplied with a ~ 1/137) are under control, the 
dominant problem at hadron colliders is the size of the logarithms (multiplied with a s ), which we link to 
the transverse momentum spectrum in the exclusive processes. 

VI. CONCLUSIONS 

Starting from charged Higgs boson production in association with a top quark, we have investigated 
processes which can be described using bottom partons. From the kinematics of the exclusive processes, 



15 



we numerically find that the factorization scale of the bottom parton has to be smaller than the threshold 
mass or the hard scale in the process: fip,b ~ M/5. In two steps we first investigate the validity of the 
asymptotic approximation in the bottom quark virtuality and then in the transverse momentum. The upper 
limit p™^ x , for which the exclusive cross section is dominated by collinear bottom quarks and large loga- 
rithms logfjOy : b/ m l)> defines the appropriate value for the factorization scale of the bottom parton in the 
inclusive process. We derive the observed dramatic decrease in the factorization scale as compared to the 
hard scale M in a process-independent approach, using only properties of the phase space and of the parton 
densities. In this simple picture we indeed find /ip jb ~ M/4. 

This observation resolves the puzzle of the discrepancy between inclusive and exclusive rates, as it 
has been present in the literature. Using an appropriate scale, the difference for example for the process 
bb — > $ [9] is not huge and well understood. Moreover, higher order calculations have been pointing to 
small bottom factorization scales [6-10]; we understand how this is caused by the partonic phase space 
without relying on numerical arguments based on higher order corrections. 

Turning around the argument, we can specify how large the logarithms actually are, which are resummed 
in the bottom parton picture. Again, they are smaller than the naive guess log(M 2 /mf) would indicate. In 
particular, for a light neutral Higgs boson one can debate using the (resummed) bottom parton cross section, 
or just integrating over the exclusive cross section [33]. For heavier neutral or charged Higgs bosons the 
logarithms are certainly large enough to require a resummation. 
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